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Abstract 

Pair production of new vector-like quarks in pp collisions is consid¬ 
ered model independent as it is usually dominated by QCD production. 
We discuss the interpretation of vector-like quark searches in the case 
that QCD is not the only relevant production mechanism for the new 
quarks. In particular we consider the effect of a new massive color octet 
vector boson with sizeable decay branching ratio into the new quarks. 
We pay special attention to the sensitivity of the Large Hadron Col¬ 
lider experiments, both in run-1 and early run-2, to differences in the 
kinematical distributions from the different production mechanisms. 
We have found that even though there can be significant differences 
in some kinematical distributions at the parton level, the differences 
are washed out at the reconstruction level. Thus, the published exper¬ 
imental results can be reinterpreted in models with heavy gluons by 
simply rescaling the production cross section. 



1 Introduction 


The discovery of the Higgs boson mm at the Large Hadron Collider (LHC) 
has given ns, for the first time, direct experimental access to one of the 
least nnderstood aspects of elementary particle interactions, the mechanism 
of electroweak symmetry breaking (EWSB). The absence of a dynamical 
mechanism that sets the scale of EWSB in the standard model (SM) makes 
it difhcnlt to conceive that the SM is the whole story. It is far more likely 
that there is new physics beyond the SM reponsible for the dynamical re¬ 
alization of EWSB, and a great deal of effort from the LHC experimental 
collaborations is going into searching for this possible new physics. The in- 
gennity of model-bnilders has provided ns with an enormons variety of new 
physics models, with a wealth of different phenomenological signatnres, to be 
searched for at the LHC. The experimental collaborations have responded 
to this theoretical richness in a three-fold way. First, they have made a 
tremendons effort in trying to cover all possible experimental signatnres of 
known models. Second, they are providing more and more details of their 
experimental analyses, thns making the recasting of snch analyses for dif¬ 
ferent models easier. Finally, they try to interpret their resnlts in terms of 
the largest possible nnmber of different models so that they are as widely 
applicable as possible. A good example of this are the searches for new 
vector-like qnarks (VLQ) that have gone from being searched for in only 
a few final states and interpreted in terms of a single model to combining 
all possible final states and reporting the corresponding bonnds in terms of 
arbitrary branching ratios (BR), see for instance |3ll3j. 

The last point mentioned above, the interpretation of experimental re¬ 
snlts in terms of a limited nnmber of theoretical models, remains one of the 
main bottlenecks in the crncial exchange of information between the theoret¬ 
ical and experimental commnnities in the search for physics beyond the SM. 
Important efforts have been made recently, like the nse of simplified mod¬ 
els [5] or the interpretation of the experimental resnlts in terms of effective 
Lagrangians. However, even this does not always captnre the relevant infor¬ 
mation for some well-motivated models, either becanse they are too complex 
to be cast in the form of a reasonably simple model or becanse the inter¬ 
pretation in terms of effective operators reqnires too many free parameters 
or is complicated by the presence of new light particles in the spectrnm. 
The recent development of Monte Carlo and fast detector simnlation tools 
has simplified in a great manner the reinterpretation of experimental resnlts 
in terms of new models. However, a serions analysis beyond an edncated 
estimate reqnires a significant amonnt of compnting power and it will never 
reach the complexity and precision of a fnll-fledged experimental analysis. 

When a new physics model has experimental signatnres close to the ones 
that are being searched for at the LHC it is worth stndying in detail how 
the corresponding analyses can be reinterpreted in terms of the new model. 
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There are essentially two possibilities, either the kinematical distributions 
relevant for the analyses are similar enough (within the experimental preci¬ 
sion) to the ones in one of the models that are being used for the interpre¬ 
tation or some of the distributions are significantly different. In the former 
case the reinterpretation is trivial whereas in the latter it becomes much 
more complicated but on the other hand there is a possible handle to dis¬ 
tinguish the new model from the ones being used for interpretation (and 
to increase the sensitivity of the search to the new model). The ideal case 
is actually somewhere in between, with identical distributions for the vari¬ 
ables that have been used in the experimental analyses but with significant 
differences in extra variables that could be eventually used to design more 
sensitive searches. The main goal of this article is to perform such a study 
in a very well motivated model. In particular we consider the reinterpreta¬ 
tion of VLQ searches in terms of a new production mechanism, namely the 
s-channel exchange of a massive color octet vector boson (a heavy gluon), in 
composite Higgs models (CHM) with partial compositeness. 

CHM [6] with partial compositeness [7] predict new VLQ, called top 
partners, that, in the light of the Higgs mass, are expected to be relatively 
light |8j, with masses in the TeV region. The mechanism of partial compos¬ 
iteness, that implies a linear mixing between the elementary and composite 
sectors, guarantees the presence of colored resonances in the spectrum and in 
particular, massive color octet vectors are naturally expected. If these heavy 
gluons have similar masses to their electroweak counterparts, they have to 
be in the multi-TeV range (due to stringent constraints from electroweak 
precision data on the latter M) and therefore it is likely that their decay 
into top partners, which is favored by large couplings, is kinematically open. 
In this case the heavy gluon will have a sizeable width and will decay mainly 
into the top partners, in pairs if kinematically allowed, thus significantly 
reducing the sensitivity of tt resonance searches to the heavy gluon jloQ 
In this case VLQ searches can be the most sensitive probes of the heavy 
gluon (see m for a recent study in which bounds on the heavy gluon have 
been computed by recasting experimental searches for pair production of new 
charge 5/3 VLQ and [T3] for implications on recent reported excesses at the 
LHC.). This is therefore an example in which we have the same final state 
considered in several experimental searches (VLQ pair production) but with 
a new production mechanism. The goal of our study is to investigate if the 
large mass of the heavy gluon has a significant impact on the kinematical 
distributions used in VLQ seaches. We will see that the relevant distribu¬ 
tions, which depending on the region of parameter space can be significantly 
different at the parton level, are very similar to the ones in which only QCD 
production is considered, once detector and reconstruction effects are taken 

^The case in which the heavy gluon decays into a SM quark and a VLQ has been 
studied in [111 112) . 
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into account. Thus, in most of the parameter space, VLQ searches can be 
trivially reinterpreted in terms of a rescaling of the cross section in models 
with heavy gluons without the need of costly computer simulations. One 
can simply use the published experimental limits on the VLQ production 
cross section as a function of the VLQ mass, that assume QCD production, 
and overlay the theoretical production cross section in the model with heavy 
gluons to compute the limits. We then use this approach to compute the 
current bounds on models with a heavy gluon as derived from the LHC run-1 
VLQ searches and estimate the sensitivity of the early LHC run-2 data. 

The rest of the article is organized as follows. We introduce the model 
in Section]^ We compare the relevant kinematical distributions in a model 
in which QCD is the only production mechanism with the ones in our model 
in Sections 3.1 at truth parton level and 3.2 including detector and recon¬ 
struction effects. We then discuss the effect that the differences between the 
kinematical distributions in both models have on the actual limits extracted 
in Section |3.3| Section is dedicated to present the current limits on the 
model and early run-2 expectations and we conclude in Section]^ Technical 
details of the model are presented in Appendix [A| 


2 The Model 


CHM with partial compositeness represent a well-motivated class of models 
of dynamical EWSB. They usually predict a number of relatively light VLQ 
with large couplings to the top (and bottom) quark, the so-called top part¬ 
ners, and also new vector resonances, possibly in the octet representation of 
the color group that couple sizeably to the top quark and its partners. They 
are therefore a prime candidate to study the interplay of VLQ and heavy 
gluon searches that we mentioned in the introduction. For the sake of con¬ 
creteness we have considered a simplihed version of the Minimal CHM m, 
based on the S'0(5)/S'0(4) coset, with a fully composite right-handed (RH) 
top quark [16]. In particular we take the MCHds model in which the SM 
left-handed (LH) quarks are embedded in a 5 of SO(5) and the top partners 
span a 4 representation of SO(4). Specihcally they read, in the basis we are 
considering. 
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where I = 1,...,5 and i = 1,...,4, respectively. In terms of SU{2)l x 
U{1)y representations, the 4 of SO(4) T gives rise to two doublets with 
hypercharges 1/6, {T,B), and 7/6, (V 5 / 3 ,A 2 / 3 ), respectively. The latter 
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contains an exotic state with charge 5/3, -^ 5 / 3 , and a charge 2/3 state, ^ 2/3 
whereas T and B have electric charge 2/3 and —1/3, respectively. After 
EWSB, a linear combination of T and ^ 2 / 3 ) that we denote ^ 2/3 
degenerate with X^j^. The orthogonal combination, that we call T', and B 
are somewhat heavier with a small mass splitting. In most of the parameter 
space their decay BR read 


RR(A5/3 ^ 

tw+) 

= BR{B ^tw-) = 1, 

(2) 

BRiX',/, 

tZ) 


(3) 

BR{T' 
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(4) 


The new ingredient with respect to the MCHds model is the presence 
of the heavy gluon. We introduce it following the partial compositeness 
mechanism m- The relevant part of the Lagrangian, in the elementary- 
composite basis, reads 


C = 


+ 


qiBXqi + tRxEftR + + i/)'!' — 

\ci{^ R)a^^d!-^tR + yfiQlYUii^R + yc2f{QlYUi5tR + h.c. 

(5) 


This Lagrangian is determined by three dimensionless (order one) couplings, 
y and ci^ 2 ) the strong couling scale /, the composite gluon mass Me and 
the ratio of elementary to composite gluon couplings Qc- The explicit 
expressions of the different operators in Eq. Q are given in Appendix 
The relevant feature for us is that the elementary gluon only couples to 
elementary quarks, with coupling and the composite gluon only couples 
to the composite quarks with coupling gc- (For simplicity we have assumed 
a universal coupling of the composite gluon to the composite states. This 
assumption does not have any relevant implication on the main conclusions 
of our study.) The elementary-composite gluon system can be brought to 
the physical basis by means of the following rotation 


(^ (cos 03 - sin 03 \ / \ 

\Gl) Vsin03 cos03yVG'J’ 


( 6 ) 


where the ratio of couplings fixes the mixing angle tan 03 = After this 
rotation we have a massless color octet, the SM gluon 5 ^, and a heavy gluon, 
G^, with mass 


Mg = 


Me 


cos 03 

The SM gluon couples universally with coupling strength 


(7) 


gs = 5 c sin 03 = 5 eCos 03 , 


( 8 ) 
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whereas the heavy gluon has couplings to elementary and composite helds 
given by 

g2 

G'^elemV’elem ■ ; ^ ^; G^V^compV^comp • \/9c 9s- (9) 

V 9c 9s 

In the following we hx the input parameters to the following values: 

g^ = 3, f = 800 GeV, ci = 0.7, C 2 = 1.7 , (10) 

where the value of gc is fixed only to test the independence of the relevant 
kinematical distributions of the reconstructed objects on the presence of the 
heavy gluon, described in detail in the next section. In order to get the 
bounds on Mq from current data and the expected reach with the early 
LHC run-2 data we will vary 2 < gc < 5. y is hxed by the top mass and 
its value determines the degree of compositeness of the LH top. It is most 
sensitive to the value of C 2 and has a milder dependence on the mass of the 
top parners. For the values of the input parameters chosen above it ranges 
from y « 0.95 for Mq, = 600 GeV to y ^ 0.61 for Mq, = 1.6 TeV. This in 
turns corresponds to a degree of compositeness for the LH top quark 



of sl = 0.78 and sl = 0.29, respectively. As it will become clear in the 
following sections, our goal is to show that the kinematical distributions of 
the reconstructed objects in VLQ searches are not affected in a signihcant 
way by the presence of the heavy gluon. This effect is not very sensitive to 
the degree of compositeness of the LH top quark, a fact that we implicitly 
check by not keeping sl fixed in our analyses. 

The only remaining parameters are the values of the top partners and 
heavy gluon masses, Mq, and Mq, and the behaviour of the light generations. 
Regarding the masses we will vary them in the relevant range 

600 GeV < Mvi, < 1.6 TeV, 1.5 TeV < Mq < 4.5 TeV. (12) 

Finally, for simplicity we assume that the hrst two generations and the RH 
bottom quark are purely elementary. 

One hnal comment to fully characterize our model is the fact that, due to 
the large number of decay channels into top partners and the large coupling 
to them, the width of the heavy gluon is typically very large, becoming quite 
easily comparable to Mq itself. In that circumstance, it is compulsory to 
consistently include the full energy-dependent quantum width HD Ea Il3|. 
Due to the rapidly decaying parton distribution functions (PDF), a larger 
width for the heavy gluon makes it more similar to the continuum QGD 
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production of the top partners than a narrower one. Thus, since we are 
trying to hnd how different the VLQ pair production via the heavy gluon is 
from the QCD one, we have considered a modihcation of the model in which 
the couplings to the composite fermions, including the RH top and the LH 
top and bottom, are rescaled down, when necessary, to hx the maximum 
width of the heavy gluon to 20% of its mass. We emphasize again that this 
variation goes in the direction of making the eventual differences between the 
two production mechanisms larger and therefore if no signihcant differences 
are found in this case, they are not expected to arise in the usual case with 
a wider heavy gluon. 


3 Comparison of kinematical distributions and rein¬ 
terpretation of limits 


In this section we want to investigate whether the presence of a heavy gluon 
in CHM affects VLQ pair production in a signihcant way. If it does, detailed 
simulations are needed to recast the result of current VLQ searches for the 
model at hand, although then new experimental searches might be designed 
to gain further sensitivity to these models. If it does not, the bounds can 
be trivially reinterpreted in CHM with heavy gluons by simply rescaling the 
corresponding cross section. In the case of a heavy gluon, one could expect 
that its large mass provides the VLQ with a larger energy and therefore 
modify the kinematical distributions of the particles they decay into. Clearly, 
how large the effect is depends on the relation between the heavy gluon and 
the VLQ masses and also on the available energy. 

In order to test all these features in detail we will concentrate on one par¬ 
ticular analysis, the VLQ pair-production searches in the Zt channel |20) and 
display the kinematical distributions assuming production only of In 

Section [3T] we will show the differences on the main kinematical distributions 
at the parton level. For a better comparison we will display the distributions 
both normalized to their corresponding cross section and also to unit area. 
This latter normalization gives us a more clear picture of whether the dif¬ 
ferences can be interpreted as purely cross section rescaling or they present 
signihcant shape differences. We will then compare the distributions after 


hadronization, showering and detector simulation is included in Section 3.2 


Finally, in Section 3.3 we compare the current limits from run-1 data with 
a full recasting of the analysis in the model versus a simple rescaling of the 
corresponding cross sections. 

We have generated a UFO model [21] using FeynRules |22| . The sim¬ 
ulations at parton level have been done with MG5 |23| while for hadroniza- 


^With this modification, the heavy gluon has a behaviour more similar to the narrower 
electroweak vector boson resonances, recently studied in m- These resonances might be 
responsible for the reported diboson excess if the decay into top partners is not open m- 
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tion/showering and detector simulation we have used pythia 6 |24) and 
Delphes 3 |25], respectively. We have used the default ATLAS Delphes 
card with the following modified parameters: the FastJet |26) ParameterR 
and the AR in the 6-jets have been set to 0.4; the 6-tagging efficiency has 
been set to 70% and the mistag probabilities for c-jets and light-quark/gluon 
jets to 20% and 0.7%, respectively. 

3.1 Kinematical differences at parton level 

In this section we explore the kinematical differences, at the parton level, 
between the production of a pair of vector-like quarks via QCD and consid¬ 
ering the presence of a heavy gluon in the s-channel. In order to disentangle 
the origin of any possible kinematical differences, we compare the pair pro¬ 
duction of quarks via QCD (QCD), via the heavy gluon only (HG) and 
via both QCD and heavy gluon, including interference effects (QCD+HG). 
We have studied a large number of kinematical distributions for different 
values of the G and masses, both at y/s = 8 and 13 TeV and we will 
report on the most relevant ones. When showing kinematical distributions in 
this section, we will show the distributions normalized to their corresponding 
cross section on the left and the same distributions normalized to unit area 
on the right. The former give us an idea of the size of the different con¬ 
tributions whereas the latter give us intuition on how the efficiencies might 
change depending on the production mechanism. 

In experimental searches it is common to start looking at distributions 
of objects in the final state that are easy to reconstruct. An example would 
be the pT of the Z boson in the Z decay of the Xy^ i^ 2/3 Zt —)• ZWb). 
We show the pt{Z) distribution in Fig. for ^/s = 8 and 13 TeV at the 
top and bottom rows, respectively. As mentioned above, the left plots are 
normalized to the corresponding cross section whereas the right ones are 
normalized to unit area. We have fixed Mq = 3.5 TeV and = 1 TeV 

to have a sizeable contribution from the HG channel. It can be seen that, 
even though there are clear differences between QCD and HG, the shapes 
of QCD and QCD+HG are very similar, both at 8 and 13 TeV. We have 
observed a similar behaviour in properties of objects near the detector such 
as the pt of leptons or other vector bosons present in the decay chain. 

As one could expect, kinematical distributions of objects further up the 
decay chain are more sensitive to the presence of the heavy gluon. We show 
in Figs. ^ and 1^ the distribution of p^iXy^) and M{Xy^X''^^^), respectively. 
The same layout is used with data at 8 (13) TeV in the top (bottom) row. The 
left panels normalized to the corresponding cross sections whereas the right 
panels are normalized to unit area and the same values of the G and Xy^ 
masses are chosen. As we see in the figures, the transverse momentum of the 
VLQ already shows some differences in shape between QCD and QCD+HG, 
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(a) 


(b) 




(c) 


(d) 


Figure 1: Distribution of the Z boson transverse momentum for -y/s = 8 
TeV (top) and 13 TeV (bottom) at parton level with Mq = 3.5 TeV and 
= 1 TeV. The plots on the left are normalized to the corresponding 
cross sections whereas the ones on the right are normalized to unit area. 


differences that, as expected, are even larger for the invariant mass of the 
VLQ pair. The differences are also larger at 13 TeV than at 8 TeV as can 
be easily understood from Fig. Due to the PDF suppression, a heavy 
gluon with Mq = 3.5 TeV is largely produced off-shell for \fs = 8 TeV. At 
13 TeV, on the other hand, the heavy gluon is mostly produced on-shell, with 
a distinctive peak in the invariant mass distribution. The kinematical distri¬ 
butions of its decay products are therefore harder and easier to distinguish 
from the QCD case. Of course, there is a balance between distributions that 
are more easily measured but less sensitive and those that are more sensitive 
to the presence of the heavy gluon but more difficult to measure. Further¬ 
more, the latter will naturally suffer a larger degradation when more realistic 
effects, like hadronization and detector simulation are taken into account, as 
discussed in Section [3.21 

Before considering detector simulation we would like to point out that 
the presence of a heavy gluon does not always imply a harder spectrum 
than in QCD production. In the plots above, we chose the Mq = 3.5 TeV 
and = 1 TeV to have a comparable contribution from QCD and 

HG productions with a VLQ mass compatible with current experimental 
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Figure 2: ^ 2/3 ^'T distribution for ^/s = 8 (top) and 13 TeV (bottom) at 
parton level with Mq = 3.5 TeV and = 1 TeV. The plots on the left 

column are normalized to the corresponding cross section whereas the ones 
on the right are normalized to unit area. 


bounds. In this example the heavy gluon mass is large enough to allow the 
VLQ to always be produced on-shell. An interesting situation occurs when 
the heavy gluon mass is slightly below the kinematical threshold for the 
VLQ pair production and the HG production dominates over the QCD one. 
In this case, there is still a sizeable production of VLQ pairs in which at 
least one of the two VLQ is slightly off-shell, thus making the corresponding 
spectrum softer than the one in QCD production. An example of this can be 
seen in Fig. |^in which we show the P't{X 2 ^^) (top) and Mx'^ V' 2/3 (bottom) 
distributions for Mq = 2.5 TeV and Mx^^^ = 1.3 TeV at = 8 TeV. 
As usual the left plots are normalized to their corresponding cross sections 
and the right ones to unit area. We see in the hgure that the heavy gluon 
resonance slightly below the pair production threshold makes the spectrum 
softer. As a result of this the corresponding analysis efficiencies are smaller 
for the QCD+HG production than the ones for QCD. This effect continues 
at ^/s = 13 TeV despite the larger available phase space. 
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Figure 3: X'^i^X' 21 ^ invariant mass distribution for y/s = 8 (top) and 13 
TeV (bottom) at parton level with Mq = 3.5 TeV and = 1 TeV. The 

plots on the left column are normalized to the corresponding cross section 
whereas the ones on the right are normalized to unit area. 

3.2 Kinematical differences after detector simulation 

Once we have explored the sensitivity of different kinematical distributions 
to the presence of a heavy gluon at the parton level, we consider now a more 
realistic set-up in which showering/hadronization and detector simulation 
are included. In order to produce the kinematical distributions after detector 
simulation the physics objects need to be properly reconstructed. We will 
describe our reconstruction method in detail in the next section and report 
here simply the corresponding kinematical distributions. Our goal is to show 
how the distributions used as example in the previous section get modified 
when detector simulation is included. However, aiming at the recostruction 
of the full X'^i^X' 2 iz system is too ambitious if one is interested in run-1 or 
early run-2 data. Thus, instead of the X'^i^X' 21 ^ invariant mass distribution, 
we will show another distribution that naturally appears at the detector level, 
namely the scalar sum of the transverse momenta of all jets in the event, //t- 
We show in Fig. the invariant mass of the objects that reconstruct the 
VLQ (top left), to show that our reconstruction method indeed captures the 
correct final state objects; the transverse momentum of the Z bosons (top 
right); the transverse momentum of the reconstructed VLQ (bottom left) and 
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Figure 4: ^ 2/3 transverse momentum (top) and Xy^X' 21 ^ invariant mass 
(bottom) distributions for ^/s = 8 TeV with Mq = 2.5 TeV and = 

1.3 TeV. The left plots are normalized to the corresponding cross sections 
whereas the right ones are normalized to unit area. 


H'y (bottom right). In these plots we have fixed Mq = 3.5 TeV, = 1 

TeV and ^/s = 13 TeV and we show only distributions normalized to unit 
area. The distributions at 8 TeV are similar to the ones we show, if anything 
with smaller differences between the QCD and QCD+HG productions. Thus, 
we see that the somewhat small differences seen at the parton level are 
almost completely washed out once detector simulation is included. This 
conclusion remains valid for all the kinematical distributions that we have 
studied. In the case that the heavy gluon is slightly below the kinematical 
threshold for pair production and dominates over the QCD production, in 
which we saw noticeable shape differences between the different production 
mechanism we have checked that these differences are still visible at the 
detector level. These kinematical differences are however overwhelmed by 
the large enhancement in the cross section and have little impact on the 
final limits. 
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Figure 5: Kinematical distributions after detector simulation for ^/s = 13 
TeV with Mq = 3.5 TeV and Mv’ = 1 TeV. The distributions shown are 

^ 2/3 

the invariant mass of the objects reconstructing ^ 2/3 fbe Z boson 

transverse momentum (top right), the reconstructed transverse momen¬ 
tum (bottom left) and H'y (bottom right). All distributions are normalized 
to unit area. 


3.3 Recasting VLQ searches 

As we have seen in the previous section, the small kinematical differences 
observed at the parton level between the QCD and QCD+HG production 
mechanisms are almost completely washed out once detector simulation is 
included. This means that the efficiencies of the experimental analyses, and 
therefore the experimental sensitivity, are likely to be very similar in mod¬ 
els with and without a heavy gluon in the spectrum. In this section we 
are going to perform a detailed quantitative analysis of the impact that the 
small remaining differences have on the final experimental limits. We will 
do so by computing the expected limits in the Mq — plane in two 

different ways. The first way is to compute the limits by performing the 
corresponding analyses on the actual model with both ^ 2/3 
second way is to compute the experimental sensitivity of the analyses as¬ 
suming only QCD production, as currently done by the LHC experimental 
collaborations. In this second approach the experimental limit on the corre- 
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Selection 

QCD cut efficiency (%) 

QCD+HG cut efficiency (%) 

Leptonic Z 

1.28 

1.36 

> 2 jets 

99.82 

99.90 

> 2 5-jets 

64.25 

64.00 

Pt{Z) > 150 GeV 

93.09 

92.50 

Ht > 600 GeV 

94.19 

93.42 


Table 1: Event selection and cnt efficiencies for both QCD and QCD+HG 
prodnctions at ^/s = 8 TeV for Mq = 3.5 TeV and = 1 TeV. Each 

efficiency is derived based on the nnmber of events which passed the cnt 
before. The efficiency of the hrst cnt inclndes the leptonic BR of the Z. 


spending cross section is a fnnction of only. We will then overlay the 

theoretical cross section in onr model for different valnes of Mq, assnming 
that the efficiencies are identical in both cases and obtaining in this way the 
limits as a fnnction of Mq- The hrst approach corresponds to the properly 
compnted limits, specihc for the model at hand. The second approach, on 
the other hand, relies on the (reasonable, as we have seen) assnmption that 
the experimental efficiencies are quite insensitive to the presence of the HG 
bnt has the advantage that it nses the experimental information that LHG 
collaborations are cnrrently pnblishing and simply rescaling the theoretical 
prodnction cross section. 

In order to fnlly exploit the small kinematical differences observed after 
detector simnlation we will apply the procednre described above to two dif¬ 
ferent experimental analyses. The hrst one is a recast of a search for pair 
prodnction of VLQ with at least one of them decaying into a Z pnblished by 
the ATLAS collaboration m, in which the VLQ is not fnlly reconstrncted. 
The second one is a more sophisticated mnltivariate analysis, to take fnll 
advantage of the small kinematical differences, in which the VLQ if fnlly 
reconstrncted. Let ns discnss the two analyses and the corresponding limits 
in tnrn. 

In onr hrst analysis we have replicated the ATLAS selection for the dilep- 
tonic channel which is snmmarized in Table together with the efficiency 
of each cnt at 8 TeV for Mq = 3.5 TeV and -^x'^g = 1 TeV. First, a pair of 
opposite-sign same-havonr leptons (pt(^) > 25 GeV, \rj{e)\ < 2.47, exclud¬ 
ing electrons within 1.37 < \r]\ < 1.52, and |r/(/r)| < 2.5) are reqnired with 
\mi£ — mz\ < 10 GeV. The events passing the hrst cnt are reqnired to have 
at least 2 jets {ptU) > 25 GeV, \rj{j)\ < 5), with at least two of them being 
6-tagged. Two fnrther cnts are applied to rednce backgronnd contamination: 
PTiZ) > 150 GeV and iLx > 600 GeV, where is the scalar snm of the pT 
of all selected jets. As we can see in the table the efficiencies are very similar 
for the QGD and QGD+HG prodnctions. We have also simnlated the two 
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Figure 6: Signal and background distributions at y/s = 8 TeV corresponding 
to (left) for events passing the first three cuts and and the invariant 
mass of the Zb system applying all selection (right). The distributions are 
normalized to unit area. We have fixed Mq = 3.5 TeV and = 1 TeV. 


main backgrounds identified by the ATLAS analysis, namely Z+jets and ti, 
using the same tools as for the signal. For the Z+jets background we have 
generated a matched sample of Z+bb and Z+cc with up to 3 additional jets 
in the matrix element. No Z+light sample has been generated since after 
the 2 6-jet selection the contribution is negligible. The ti sample has been 
generated at NLO in MG5 and the dileptonic decay has been performed with 
Madspin |27| . 

As expected the relevant kinematical distributions after cuts present min¬ 
imal differences between the QCD and QCD+HG production mechanisms. 
As an example we show in Fig. [^the FFt distribution after the three first 
cuts (left) and the invariant mass of the Zb system (right), which is the 
discriminant variable defined as in the ATLAS publication: the Z boson 
is the candidate selected before and the 6-jet is the highest-pT 6-tagged jet 
in each event. We also show the corresponding background distibutions for 
comparison. All distributions are normalized to unit area and we have fixed 
Mg = 3.5 TeV and = 1 TeV. Finally, we show the corresponding 

limits in the Mq — Mx> plane computed with the two methods described 

2/3 

above in the left panel of Fig. The solid purple line corresponds to the 
bounds computed within the model with the HG included. The green line, 
labeled “Scaled QGD limit”, corresponds to the bounds computed assuming 
that the experimental sensitivity is the same in the QGD and QGD+HG pro¬ 
duction mechanisms. This second method is illustrated in the right panel 
of Fig. in which we show in solid black the experimental limit on the 
X'^I^X' 21 ^ production cross section as a function of assuming QGD 

production, as well as the theoretical pair production cross section, assuming 
QGD+HQ production, for different values of Mq- The points in which these 
curves cross the experimental limit determine the solid green curve on the 
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(a) 


(b) 


Figure 7: Left: 95% confidence level (C.L.) lower limits in the Mq — 
plane computed with the two methods discussed in the text. We show for 
comparison the limit reported by ATLAS {^/s = 8 TeV) in dashed red and 
the corresponding number that we have obtained with our recast of that 
analysis in dashed blue. Right: 95% C.L. upper limits on the pair production 
cross section as a function of assuming QCD production (solid black) 

with 1- and 2-cr bands in blue, together with the theoretical production cross 
sections in the model for different values of Mq- 


left panel. As we can see in the figure, the two ways of computing the limits 
give virtually identical results. This indicates that, at least with current 
analyses at 8 TeV that do not fully reconstruct the VLQ, one can simply 
use the experimental limits on the production cross section as a function of 
the VLQ mass published by the experimental collaborations to put limits on 
models with a heavy gluon. We also show in the left panel with horizontal 
lines the limit we obtain in the case of QCD production and the limit re¬ 
ported for that case by ATLAS. The level of agreement between the two and 
of them with the large Mq limit shows the consistency of our results. 

Let us now consider the second, more sophisticated, type of analysis to 
try to exploit the small kinematical differences observed at the detector level. 
In this analysis we modify slightly the event selection and reconstruct the 
^ 2^3 decay to a Z boson as follows. The Z boson is reconstructed as in the 
previous analysis. We then require at least 4 jets with at least 2 of them being 
6-tagged jets. Two of the selected jets are required to have an invariant mass 
inside a 10 GeV window of W mass and the two leptons reconstructing the 
Z boson are required to satisfy Ai?(£i,£ 2 ) < 1-5. These two last cuts reduce 
the background contamination while still leaving the px and variables 
with discriminant power. For the reconstruction of the VF, top and we 
use a to find the best combination of objects reconstructing the desired 
particles. Specifically, for a particle X which decays into two particles i,j, 
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we define Xx 


Xx = 


{M,j - Mxf 


a 




(13) 




where Mij is the invariant mass of the i,j system, Mx is the mass of the 
mother particle, ax is the width of the mother particle and AR{i,j) is the 
angular separation between both decay products. Using this dehnition we 
select the jets and 6-jet pairing that minimizes 


{xw + Xt) Ai?(Z, t). 


(14) 


Only hadronic decays of the W boson are considered and the top quark is 
reconstructed with the reconstructed W boson and a 6-jet. The AR term 
in the Xx dehnition is used to include the angular separation between the 
decay products in the object reconstruction. 

After this selection is applied we dehne a multivariate analysis using a 
neural network (NN) implemented through TMVA |28) with 7 input distribu¬ 
tions: Jet multiplicity, PT{t), pT{Xy^), pT'{Z), 

is a variable that is sensitive to both the energy of the decay products of the 
vector-like quark and the angular separation between them. It increases as 
we move to a more boosted regime and is dehned as 


aliqn 

Pt 


V^Pt(^)^ + PT{Wy + PTiby 


max (^Ai?(X'/3, Z), Ai?(A'/3, W), Ai?(X'/3,6 


(15) 


These variables are selected to provide good discrimination between signal 
and background and for their potential to provide discrimination between 
both production mechanisms (QCD and QCD+HG). The NN is trained for 
each pair of masses of G and quark separately and the output of the 
NN is used as discriminant variable to derive exclusion limits. We show in 
Fig. I the signal and background distributions of (left) and the NN 

output (right) for ^/s = 8 TeV, Me = 3.5 TeV and Mx^^^ = 1 TeV. In both 
cases the distributions are normalized to unit area. 

We use the CLg method (29) with the NN output as the discriminat¬ 
ing variable to compute the limits in the Mg — plane with the two 

methods described at the beginning of this section. The result is shown in 
Fig.g using the same color coding as in the left panel of Fig. We see from 
the results in the hgure that, although the multivariate analysis seems to be 
more sensitive to the presence of the heavy gluon, the difference between the 
properly computed limit, taking into account the presence of HG, and the 
one in which the QGD experimental limit is used is only marginal. The im¬ 
provement in the actual limits is also minimal. Thus, the conclusion drawn 
from the previous analysis that one can simply take the experimentally pub¬ 
lished limits that assume QGD production and use them to compute the 
bounds on models with a heavy gluon seems to be robust and apply also to 
more sophisticated analyses with a richer object reconstruction. 
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Figure 8: (left) and NN output (right) distributions for signal and 

background normalized to unit area for ^/s = 8 TeV, Mq = 3.5 TeV and 
Mv' = 1 TeV. 

^ 2/3 



Mo [GeV] 
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Figure 9: 95% C.L. lower limit on as a function of Mq for ^/s = 8 TeV 

with the full run-1 integrated luminosity. The green line shows the lower 
mass limit assuming that the only effect is an increase in the cross section 
and the purple line shows the lower mass limit when considering the correct 
cross section upper limit for each pair of masses. 


4 Limits on heavy gluons from vector-like quark 
searches 

Once we have checked in detail that the efficiencies and therefore the ex¬ 
perimental limits on the cross section of pair produced VLQ are insensitive 
to the presence of the heavy gluon we proceed to compute the limits in the 
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Mq — Mq plane for our CHM, where now Q stands for any of the VLQ 
present in the spectrum. In Section |4.1| we use the published run-1 limits 
on VLQ pair production (QCD only) to compute the current bounds on the 
model assuming the experimental efficiencies to be insensitive to the presence 
of the heavy gluon. We then proceed in Section 4.2 to estimate the bounds 
that can be obtained from the multivariate analysis described in the previous 
section with the early run-2 data (we assume an integrated luminosity of 10 

fb-Q. 


4.1 Current limits 

In this section we reinterpret the current results published by the ATLAS 
and CMS collaborations assuming that the experimental efficiencies are in¬ 
sensitive to the presence of the heavy gluon. We have tried to be as complete 
as possible using different analyses that target all the VLQ present in the 
model under study. 

When trying to use the published limits we find a first obstacle. We 
have checked in the previous sections the independence of the efficiencies to 
the presence of the heavy gluon, however they can be quite sensitive to the 
decay BR assumed for the VLQ. The reason is cross-contamination between 
different decay channels. A specific analysis is more sensitive to certain 
decay channels than others and therefore if the BR change we change the 
population of samples with different efficiencies and thus the global efficiency 
of the analysis. We have used a specific analysis whenever the BR assumed 
for the channel that that particular analysis is most sensitive to is equal or 
smaller than the one in our model. This gives us a conservative bound since 
we would be assuming in this way a global efficiency that is equal or smaller 
than the one in our model, due to the contamination of channels with smaller 
efficiencies. The experimental collaborations have made a tremendous effort 
in obtaining lower limits on the mass of the VLQ for arbitrary BR. In doing 
that, they have computed the corresponding limits in the production cross 
section as a function of the VLQ mass. It would be extremelly useful to 
provide that information too, which would allow us to compute the limits in 
models with heavy gluons using the correct BR. 

A second difficulty arises from the fact that the experimental limits are 
always reported assuming that only one VLQ is present in the spectrum. In 
our case, due to the specific BR pattern, we can group the VLQ in two differ¬ 
ent groups, Ai5/3 and B both decay into Wt and, unless full reconstruction of 
one decay leg with same/opposite sign leptons is performed, they are difficult 
to distinguish (see m for studies on how to characterize the different top 
partners in a similar model, assuming QCD production). A second group 
is formed by and T', which have the exact same decay patterns. In 

our model, V5/3 and are degenerate and lighter than T' and B. In the 
region of parameter space that we are exploring the mass splitting is of the 
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Figure 10: Left: Mti mass as a function of mass for the chosen 

parameters. Right: comparison between the 95% C.L. experimental limit 
obtained assuming only production (black solid line) and assuming 

production of both X^^ and T' (black dots). The solid red line represents 
the cross section of pair production of ^2/3 while the dashed green line 
represent the cross section of pair production of both ^2/3 
figure Mg = 3.5 TeV and ^/s = 13 TeV have been chosen as an illustrative 
example. Similar behaviour is found for different masses and center-of-mass 
energy. 


order of 10%, which is large enough to make interference effects small but not 
so large that two distinct peaks would show up in the relevant distributions. 
Furthermore, the heavier ones, T’ and B, are heavy enough to make their 
contribution to the total cross section much smaller than the one from the 
lighter ones, X5/3 and Xy^, in the searches to which both contribute. We 
show in the left panel of Fig. 10 the T' mass as a function of the ^2/3 ™^ss, 
for the chosen parameters. Thus, the global efficiencies are likely to be quite 
insensitive to their presence. In order to be quantitative, we have repeated 
our multivariate analysis including the two charge 2/3 quarks, X'^i^ and T/ 
for a few selected points. We show in the right panel of Fig. [^a comparison 
between the experimental sensitivity (limit on the production cross section) 
obtained in the full model with both quarks and the one obtained neglecting 
the presence of T'. As we can see the two curves are compatible within un¬ 
certainties and we can therefore compute the limits overlaying the total cross 
section (including X^^ and T') over the published experimental bounds on 
the production cross section. 

We have used the following experimental results to compute the current 
limits on the VLQ present in out model. For B we have used the results 
published in |31| (ATLAS) and |32j (CMS). The ATLAS analysis considers a 
singlet model with different BR than in our case. However, it is very sensitive 
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(a) 


(b) 


Figure 11: 95% C.L. lower limits on the Mq — (left) and Mq — 

(right) planes, derived from full run-1 published data (y^ = 8 TeV), for 
different values of Qc- The filled orange and green region represent the limits 
for CMS and ATLAS, respectively, assuming QCD production of and 
T' (left) and X5/3 (right). 


to the W decay and, as discussed, we have conservatively rescaled their result 
to our 100% BR into Wt. The CMS result is already published for a BR of 
100% to the W channel. Unfortunately, the mass of the B quark in our model 
is, for the range of parameters we have considered, typically heavier than the 
masses these analyses are sensitive to and we cannot therefore report bounds 
for all the Mq values. We can nevertheless provide the corresponding limits 
in some cases. For instance, using the ATLAS results we can exclude values 
oi Mb ranging from Mb — 750 GeV for Mq = 4 TeV to Mb — I TeV for 
Mq = 2 TeV. Using the CMS results we can excluded values ot Mb ranging 
from Mb — 800 GeV for Mq = 3 TeV to ^ 950 TeV for Mq = 2 TeV. 
These limits have been derived assuming Qc = 3. For the and T’ quarks 
we have used the results in [3] (ATLAS) and |1] (CMS). The ATLAS results 
contain a combination of two analyses that target decays of the VLQ into 
H and LF, respectively. Since charged current decays are not present in 
our model we use the result quoted for the doublet (same BR as in our 
model) only in the H channel. The CMS analysis reports limits for a 2:1:1 
BR pattern into W, H and Z, respectively, and the analysis used is very 
sensitive to the Z channels. Thus, we conservatively scale their results to 
out BR as discussed above. The results can be found in the left panel of 
Fig. 11 for different values of Qc- The corresponding value of Mb' can be 

In the case of the V5/3 the same ATLAS analysis 


obtained from Fig. 10 


used for the B quark published results for the search of a vector-like quark 
with a electric charge of 5/3 assuming a BR of 100% to Wt. We have also 
used the CMS analysis |33) searching for VLQ with an electric charge of 5/3 
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Figure 12: 95% C.L. lower limit on the ^2/3 as a function of Mq at 

y/s = 13 TeV and an integrated luminosity of 10 fb“^. The green line shows 
the lower mass limit assuming that the only effect is an increase in the cross 
section and the purple line shows the lower mass limit when considering the 
correct cross section upper limit for each pair of masses. The dashed blue 
line shows the lower mass limit when only QCD production is assumed. 


in the same-sign dilepton hnal state. They assume a 100% BR into Wt so 
we can directly use their results, which can be found in the right panel of 


Fig. 11 


4.2 Early run-2 expectations 

Once we have computed the current bounds in the Mq — Mq plane with Q 
running over all the VLQ present in our model we would like to compare 
them with the reach of a run-2 early data set. The goal is two-fold, hrst 
we would like to see how much the limits can improve with a few fb“^ of 
integrated luminosity at y/s = 13 TeV with respect to the current ones; 
second we would like to explore the sensitivity to the presence of the heavy 
gluon at the initial stages of a new run with higher energy, as opposed to 
the analyses we have performed so far in which the limits are at the verge of 
the kinematical reach. 

We take as a benchmark analysis the multivariate one we presented in 
Section 


3.2 and estimate its reach in the Mq — plane at 

and an integrated luminosity of 10 fb~^ 


s = 13 TeV 
We use the same selection and 


tools and compute the limits again in two ways, either including QCD+HG 
from the start or computing the upper limits in the cross section assuming 
QCD production and rescaling the corresponding cross section. Our results 


are summarized in Fig. 12 in which we show the 95% C.L. lower bounds 
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in mass for the -^2/3 h^ark as a function of the G mass. Two points are 
worth emphasizing. The first point is that the limits computed by a simple 
rescaling of the cross section (using the results that assume QCD produc¬ 
tion as would be published by the experimental collaborations) are perfectly 
compatible with the full ones computed with the model at hand. Thus, even 
when exploring regions of parameter space far from the kinematical reach of 
the experiment, it seems like the experimental limits on the pair production 
cross section of VLQ is quite insensitive to the presence of the heavy gluon. 
It is therefore perfectly legitimate to use them to put bounds on any model 
in which heavy gluons contribute to the VLQ production. A second impor¬ 
tant point is that there is, already with 10 fb“^ of integrated luminosity, a 
signihcant improvement in the limits with respect to current ones. The early 
run-2 limits on the mass of the VLQ are in the ~ 820 — 1160 GeV range 
(to be compared with the equivalent ~ 660 — 1000 GeV at ^/s = 8 TeV). A 
similar improvement can be expected with respect to the complete current 
limits presented in Fig. 11 Given these two points, it is clear that as soon 


as the experimental collaborations publish their results on VLQ production 
during run-2, they can be easily translated to bounds on the parameter space 
of models with heavy vectors in the spectrum. 


5 Conclusions 

We have considered in this article the problem of reinterpretation of exper¬ 
imental results in terms of models different from the ones used by the ex¬ 
perimental collaborations. Such reinterpretations are always delicate since 
it is difficult to assess the impact of the new features of the model in the 
hnal results. Assuming that the signatures of the two models are reason¬ 
ably similar, there are, broadly speaking, two options. The hrst option is 
that the distribution for the kinematical observables used in the analysis 
are very similar in both models. In that case the experimental sensitivity 
(limit on the production cross section as a function of the relevant scale) 
will be very similar for both models and the limits can be computed by a 
simple rescaling of the theoretical cross section. The alternative option is 
that some observables show signihcant differences between the two models. 
The reinterpretation of the limits in the new model become then much more 
complicated, and the only possibility is to recast the corresponding analysis 
using fast detector simulation tools. This latter case is not completely sat¬ 
isfactory, as it does not take full advantage of the experimental publication 
and requires a certain degree of mastery in the implementation of the model 
and the use of Monte Garlo and detector simulation tools. On the other 
hand, it offers extra handles on how to more efficiently search for the new 
model. 

For the sake of concreteness we have considered pair production of new 
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VLQ in CHM. A common occurence in these models is the presence of a 
heavy color octet vector boson, heavy gluon, that can contribute to the pair 
production of the VLQ. We have performed a detailed quantitative study of 
the effect that the heavy gluon has on the distributions of different kinemat- 
ical variables, both at parton level and after detector simulation. The small 
differences found at the parton level in variables that are relatively up in 
the decay chain are almost completely washed out after detector simulation. 
Thus, in practice, the limits on the VLQ and heavy gluon masses that one 
obtains by using the full model and the ones that use just QCD production 
(neglecting the heavy gluon) to compute the limit on the production cross 
section are virtually identical. The latter has the precious advantage that 
it uses directly the published limits computed by the experimental collabo¬ 
rations and require only a minimal amount of Monte Carlo simulations to 
compute the corresponding theoretical cross sections at the parton level. 

We have tested the validity of this approximation in two different regimes 
of an experiment, one in which the limits are close to the kinematical reach of 
the experiment and another in which a new energy frontier is being explored 
for the first time. These correspond to the full LHC run-1 data and the early 
run-2 data, respectively. In both cases we have found that the experimental 
limits on the production cross section are quite insensitive to the presence of 
the heavy gluon and therefore the experimental results can easily be reinter¬ 
preted in models with new VLQ and heavy gluons. This extends not only 
to current analyses but also to more sophisticated analyses in which full 
kinematical information, including the reconstruction of the VLQ, is used. 
Presumably only a full reconstruction of the VLQ pair invariant mass could 
give a significant sensitivity to the presence of the heavy gluon. 

We have reached these conclusions by performing an exhaustive set of 
comparisons of kinematical distributions and detailed statistical analyses, 
which include multivariate techniques. We have found that the main degra¬ 
dation of the HG signal occurs at the reconstruction level, due in part to 
a wrong assignment of physical objects by the reconstruction algorithms in 
the presence of boosted topologies due to the overlap of the considered jets. 
The use of boosted techniques might partially compensate this effect but a 
detailed analysis, which is beyond the scope of this study, would be needed 
to assess in a quantitative way the discriminating power of a reconstruction 
method using larger jets and jet substructure variables. In any case, our 
modeling of the detector is based on fast detector simulation and therefore 
our conclusions are limited by the inherent accuracy of such a tool. In that 
sense it would be very important that experimental collaborations consider 
the possibility of comparing the experimental bounds on the VLQ pair pro¬ 
duction cross section with and without a heavy gluon. 

We would also like to encourage the LHC experimental collaborations to 
report their limits on the VLQ pair-production cross section as a function 
of the VLQ HR. They are available to the collaborations, as they are used 
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to compute the limits on the VLQ masses for arbitrary BR. However they 
are not publicly available, thus limiting the reinterpretation of their limits 
to other models, e.g. with or without heavy gluons. 
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A Technical details of the model 


In this appendix we detail the different terms appearing in the Lagrangian 
Eq. (|^. The covariant derivatives read explicitly 

-|- g—WJj^ + + geG^f^ qi, 


— ( idfj^ + + gcCl ) tn, 

iTJn'k = ( idn H— + gcG'l^ ) 'k. 


(16) 


The remaining terms read 




{QlYUn^R 

{OiYui^tR 


^Sh[{X,/s)R}^+ - BR}¥-]tR 
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w 
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(17) 


+ 'Sfe/2^ + C^/2^2/3 


= biBR + ti 

1 


Yi/2^R + Yk/2{^2/YR 


-^^ShtltR, 


+ h.c.j 

(18) 

(19) 

( 20 ) 
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where we have denoted 


X X 

Sa; = Sm-, Cx^COSj, 


( 21 ) 


except for sw and cw, which are the sine and cosine of the Weinberg angle. 
p is the physical Higgs boson and h reads, in the unitary gauge 


h={h)+ p, 


( 22 ) 


with 

fs(^h) =v^2A6 GeV. (23) 

These terms hx the electroweak couplings of the top and its partners and also 
the top mass and mixing. In particular the latter two have some implications 
on the degree of compositeness of the top quark and therefore on its couplings 
to the heavy gluon. 
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